Astronomy & Astrophysics manuscript no. ms 


© ESO 2010 


May 3, 2010 





A diagnosis on torque reversals in 4U 1626—67 

Z. Zhangi'2 and X.-D. Li^'^ 
^ Department of Astronomy, Nanjing University, Nanjing 210093, China 

^ Key Laboratory of Modern Astronomy and Astropliysics (Nanjing University), Ministry of Education, China 
e-mail: zhangzhenl9841209(9gniail . com, lixd@nju.edu.cn 



o 

(N 

<: 
o 

6 ' 



> 

o 

in 
vn 

^' 

o 
o 



X 



ABSTRACT 

Several X-ray pulsars have been observed to experience torque reversals, which provide important observational clues to 
the interaction between the neutron star magnetic field and the accretion disk. We review the current models proposed 
for the torque reversals and discuss their viability based on the observations of the quasi-periodic oscillations (QPOs) 
in 4U 1626—67. Most of these models seem to be incompatible with the evolution of the QPO frequencies if they are 
interpreted in terms of the beat frequency model. We suggest that winds or outflows from the neutron star and the 
accretion disk may play an important role in accounting for the spin-down in disk-fed neutron stars. 
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1. Introduction 

X-ray pulsars are magnetic neutron stars (NSs) in binary 
systems, accreting from a normal companion star. They are 
excellent samples to study the mass and angular momen- 
tum (AM) tra nsfer be t ween an NS and the surrounding 
accretion flow (jNagasel |1989[ ). Mass accretion onto an NS 
may occur through either an accretion disk fed by Roche- 
lobe overflow, or capture of the stellar wind from the com- 
panion. There has been long-term monitoring of the spin 
evolution in X-ray pulsars, which is generally thought to re- 
sult from the interaction between the NS and the accretion 
flow. Especially, CGi?0/BATSE has provided continuous 
monitoring of X-ray pulsars (Bildsten ct al., 1997), among 
which 4U 1626—67 and GX 1-1-4 showed steady spin-up un- 
til a sudden torque reversal appeared with the spin-down 
rate similar to that in the spin-up stage but with oppo- 
site sign, while Cen X -3, OAO 1657-415, and Her X-1 
(|Klochkov et al.l . |2009[ ) exhibited a secular spin-up trend 
with short-time spin variations. 

Obviously the torque reversals reflect a dramatic change 
' in the pattern of interaction between the NS magnetic 
I flelds and the accretion disks. There have been lots of dis- 
1 cussions in the literature on the explanation of these in- 
triguing phenomena. However, a widely accepted model 
has not emerged. In this work we focus on the torque 
reversals in 4U 1626—67, a 7.66 s X-ray pulsar. Its op- 
tical counterpart was identified as KZ Tra, a faint blue 
star wit h little or no reddening (jMcClintock et aP . Il977t 
iBradt fc McClcntock. ig gj). T he short (42 minutes) orbital 
period ()Middleditcht al.l . Tl98lll indicates that 4U 1626-67 
is in an ultra-co mpact binary with a h ydrogen-depleted low- 
mass secondary ([Nelson et al.l . |1986[) . Observations during 
1977-1989 showed that 4U162 6-67 was spinning up at a 
rate of ~ 8.54 x lO^^^ Hz s~ MChakrabartv et all (|l997aD 
found that, from April 1991 to June 1996, the NS spun 
down at a rate of ;> ~ —7.17 x 10~^^ Hzs~^. Though 



the torque reversal was not observed directly, it was es- 
timated to occur in June 1990. More recent observations 
with Fermi/GBM and Swift/BAT showed that 4U 1626-67 
experienced a second t orque reversal near February 2008 
(jCamero-Arranz et al.l . l2010[ ). Since then it has been fol 
lowing a steady spin-up at a mean rate of ;> ~ 5 x 10 
Hzs^i. 
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While the spin changes in X-ray pulsars provide one way 
to study NS magnetic field-accretion disk interaction, quasi- 
periodic oscillations (QPOs) offer another clu e. QPOs have 
been detected in eight accreting X-ray pulsars ()Bozzo et al.l 
l2009l . and references therein), and are usually regarded as 
the signature of the existence of an accretion disk. The mea- 
sured QPO frequencies are in the range of ~ 1 — 200 mHz, 
consistent with being relevant to the inner radius of the ac- 
cretion disk around a highly magnetized NS in its bright X- 
ray state. In the most favored beat fr e que ncy model (BFM, 
lAlnar & Shaha^ . 119851 : iLamb et al.L [1985), the QPO fre- 
quency results from the beat between the orbital frequency 
at the inner edge of the disk and the spin frequency of the 
NS. Thus the QPO frequency may provide helpful infor- 
mation about the accretion rate in the disk. In this paper, 
after reviewing the current models of torque reversals in sec- 
tion 2, we comment on the feasibility of these models based 
on the BFM interpretation of the QPOs in 4U 1626-67 
(section 3). In section 4 we present a new picture for the 
explanation of the torque reversals, invoking winds or out- 
fiows from the NS and the accretion disk. We summarize 
our results and discuss their possible implications in section 
5. 



2 



Z. Zhang and X.-D. Li: A diagnosis on torque reversals in 4U 1626—67 



2. A review of proposed interpretations for torque 
reversals 

2.1. Magnetically threaded disk model 

The magnetically threa ded disk m odel was first proposed 
and inve stigated bv|Ghosh fc Lamb (1979a, b, hereafter GL. 
See also iBozzo et al.l (|2009t) for a detailed description and 
discussion about the model). Briefly speaking, in addition 
to the material torque 



No = A/(GM^o)^/^ 



(1) 



where M is the NS mass, M the accretion rate, and tq the 
inner radius of the disk, GL introduced a magnetic torque. 



N„ 



B^Sy dr, 



(2) 



where i?z and are the vertical and toroidal components 
of the magnetic field at the surfaces of the disk, respec- 
tively, Tout is the outer radius of the disk. The magnetic 
torque is produced by the twist of the magnetic field lines 
threading the disk. Thus, the pattern how the magnetic 
lines are twisted decides the scale of the magnetic torque. 
Assuming that the shear amplification of B^ occurs on a 
similar timescale to the reconnectio n of the slipp ed field 
lines, GL and late investigations (e.g. IWangj . I1995D showed 
that, 



oc 7(rJs - fix) 



(3) 



where 7 is the azimuthal pitch of the field lines, Jig and 
r^K are the angular velocity of the NS and the Keplerian 
angular velocity in the accretion disk, respectively. Defining 
the corotation radius as the radius where the Keplerian 
angular velocity of the material in the disk is the same as 
that of the NS, GL showed that the inner disk with r < Tc 
spins the NS up, since B^B^ < and A'^mag > 0, while the 
part outside Tc draws the NS to spin down. The total torque 
exerted on the NS is the sum of the magnetic torque and 
the accretion torque, which can be described as. 



(4) 



where n(ujs) is a dimensionless function of the fastness pa- 
rameter w,, 



(5) 



where B is the dipolar field strength on the surface of the 
NS. The value of n{ojs) can be either positive or negative, 
and there is a critical value of Wc, at which n(wc) — 0, 
corresponding to an equilibrium state. If the NS is a slow 
rotator, lOg < Wc, then n(a;s) > 0, and the NS experiences a 
spin-up torque; when Wc < Wg < 1, n(ws) becomes negative. 

In the GL-type model, if M does not change much, given 
sufficiently long time of evolution, the NS will always reach 
the equilibrium state. The most straightforward explana- 
tion for the torque reversals in X-ray pulsars is the change 
of M . However, this explanation has met several difficul- 
ties: (1) for Cen X-3, a fine-tuned change of M is required 
(jBildsten et al.L 119971 ). while observations suggested that 
the mass accretion ra te did not vary consid erably during 
the torque reversals (|Raichur fc Paul I2008D : (2) for GX 



1+4, the X-ray lu minosity seemed to increas e during the 
spin-down episode ("Chakraba rtv et all Il997bf ) , which is in 
contrast with the prediction of the original GL model. 

In a modified version of the magnetically threaded disk 
model. iTorkelssonI ()1998D suggested that, in addition to the 
GL-type torque, there is a torque coming from the inter- 
action between the stellar magnetic field and the disk's 
own, dynamo-generated magnetic field. The azimuthal field 
strength B^ is shown to be cx r"^/"*, decreasing with in- 
creasing r more slowly than the formalism i?^ oc sup- 
posed by GL. Therefore, the latter form of the magnetic 
torque dominate at large radii. The total torque is then, 

A^tot = A^GL ± A^mag,dyn. (6) 

Whether the dynamo-generated magnetic torque is positive 
or negative is not decidable, since the direction of the disk 
field is arbitrary, and is not determined by the difference in 
angular velocity between the star and the disk as in GL, but 
a transition between them is reflected as the torque reversal, 
if the magnetic field in the disk reverses. Without invoking 
definite relation between M and the spin state, Torkelsson 
(1998) indicated that the timescale for an NS staying in one 
stable spin state is uncertain, but the timescale for revers- 
ing the torque is the same as that for reversing the magnetic 
field via diffusion, or the viscous timescale. This might ex- 
plain the short timescale for the reversals compared with 
the long timescale over which the torque remains the same 
as the difference between the diffusive (viscous) timescales 
for the inner region of the disk and the entire accretion 
disk. The complication of the model lies in numerous pa- 
rameters, among which quite a lot are unclear, so it is not 
certain how important the dynamo-generated part is in the 
whole expression. Besides, this model fails to explain the 
stability of the torque (over a timescale of several years) in 
4U 1626-67. 



2.2. Retrograde or warping disk model 

If a prograde disk provides the NS with positive torque, 
why cannot a retrograde disk spins the NS down? In the 
simplest situation, the torque the star receives comes from 
the accreted material, so that 



A^ = ±M(G'Mro)^/^ 



(7) 



corresponding to a progarde and retrograde disk, respec- 
tively. A transition between prograde and retrograde rota- 
tion of the disks presents as the torque reversal naturally. 
Moreover, from Eq. (7), a similar spin- up and spin-down 
rate is anticipated, just coincident with the observations of 
the spin evolutions of GX 1-^4 and 4U 1626-67. 

iNelson et al.l (jl997^ re freshed the idea of the retrograde 
disk for GX 1-1-4 (Makishi ma et al.l . [19881 ). considering the 



fact that the X-ray luminosity increases during the spin- 
down episode. As described above, with the help of the ret- 
rograde disk, the X-ray luminosity can increase both in the 
spin-up and spin-down episodes, which successfully avoids 
the problem faced by the GL-type model. 

While ret rograde disks may ex i st in wind-fed systems 
like GX 1-H4 (iMats uda et al.lll987Lll99ll:jB6rner fc Anzerl 
[T99I lRuffertl .[T999: KrvukoveLaD, |2Q05|) , for Roche-lobe 
overflow systems like 4U1626— 67, the specific AM initially 
carried by the accretion stream is comparable to the spe- 
cific orbital AM of the companion star, and should circu- 
larize in the prograde sense well before reaching the NS 
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magnetosphere (Lubow & Shu 1975). This problem might 
be solved by disk warping induce d by X-ray irradiation 
from the central sourc e (Pringle, Il996f) or by the mag- 
netic torque (|Lail . I1999D . If the inner part of the accretion 
disk flips over by more than 90°, and rotate in the oppo- 
site d irection ()van Kerkwiik et all . 119981 : IWiiers fc Pringlel 
Il999t ). this would lead to a torque reversal. However, it 
is highly uncertain whether irradiation can cause such a 
warping and whether the warped disk can remain for a suf- 
ficiently long time. More recent theoretical work and three- 
dimensional simulations showed only slight 10 — 20%) 
warping of the disk around misalig ned magnetic stars 
(|Terauem fc Paoaloizou . .2000: .Romanova et al.L l2003[) . 



2.3. Propeller model 

A general requirement for stable accretion from a disk to 
an NS is that the velocity of the NS magnetosphere must 
be smaller than the local Keplerian velocity at the inner 
edge of the disk, otherwise the propeller mechanism will 
prohibit f urther accretion by ejecting the accreted mate- 
rial away (jlllarionov fc Sunvaevl . [l975[ ). The propeller mo- 
tion occurs at the boundary of the magnetosphere, so the 
ejected material carries the AM away from the NS, and 
spins the NS down. One problem related to the propeller 
model for the torque reversals in 4U1626— 67 might be 
that X-ray pulsations were detected during the spin-down 
stage, indicating that accretion was still going on eve n when 
ro > Tc, though in simulat i ons b y lUstvugova et al.l (|2006[ ) 
and I Romanova et"al] ()2004 [20051 ) both accreti on and spin- 
down were observed at the propeller stage. iPerna et al.l 
(|2006[ ) suggested a model for simultaneous accretion and 
ejection around magnetized NSs. When the spin axis of an 
NS is not aligned with the magnetic axis, the inner radius 
of the disk relies on the tilt angle between the two axes 
and the longitude 0, so that the inner edge of the disk is 
not circular. For certain values of 9, the inner disk radius 
is partially larger than where the propeller mechanism 
starts up. Meanwhile, the NS accretes in the region where 
the inner disk radius is smaller than Tc- In other words, the 
system can undergo the propeller and accretion phases at 
the same time. Moreover, a fraction of the ejected material 
does not receive enough energy to be completely unbind, 
and hence falls back into the disk. Thus it is possible that 
for a given accretion rate of the NS, there are multiple so- 
lutions of the mass flow rate through the disk. The spin 
evolution is determined by the AM transferred from the 
disk to the NS through accretion, and that given by the NS 
to the ejected matter. When is larger than a critical value, 
the system may settle in a limit cycle of spin- up/spin-down 
transitions for a constant value of the mass accretion rate. 



2.4. Bi-state model 



(ADAF) with small changes in the accretion rate. When 
M becomes smaller than a critical value Mcr, the inner 
part of the accretion disk may make a transition from a 
Keplerian, thin disk to a sub-Keplerian ADAF, in which 
the angular velocity in the disk r2(r) = Ail,K{r) with A < 1 
(iNaravan fc Yj FT995h . In this case the corotation radius 
becomes r[. = A^I^Tc < Tc (Here we use the prime to de- 
note quantities in ADAF). The torques exerted on the NS 
by a Keplerian and sub-Keplerian disk can be estimated to 
be 



l-(8/7)(ro/re)^/^ 
e " l-(ro/rc)3/2 ' 

and 

7^,l-(8/7)(r^,/r^)3/2 



1 - {ryr^y^ 



(8) 



(9) 



respectively, where Nq = A/(G'A/rQ)^/^. Since the inner ra- 



dius of the disk does not change much (ro 



'01' 



the dy- 



namical changes in the disk structure may lead to the slow 
{ujs < i^c, spin- up) and rapid (cjg > Wc, spin-down) rotator 
sta ge alternatively when M varies around Ma. 

iLovelace et al.l ()1999D developed a model for magnetic, 
propeller- driven outflows that cause a rapidly rotating mag- 
netized NS accreting from a disk to spin-down. An impor- 
tant feature of their results is that the effective Alfven ra- 
dius TA depends not only on M and B, but also on fig. 
Because ta decreases as ilg decreases, there exists a mini- 
mum value of fig for stable accretion disks, and for a given 
f2s, there could be two values of ta, one larger than Tc and 
the other smaller than re- This points to a mechanism for 
the propeller from being "on" to being "off' , when there is 
a change between the two possible equilibrium configura- 
tions, leading to transitions between spin-down and spin- 
up with roughly similar rates for nearly constant M. Since 
the transitions may be stochastic, and triggered by small 
variations in the accretion flow or in the magnetic fle ld con- 
figuration, this model, similar as iTorkelssonI (|1998[ ). could 
be responsible for the torque reversals in Cen X-3 rather 
4U 1626—67, provide d th at ta ^ '"co- 

iLocsei fc Melatod (|2004l) presented a disk- 
magnetosphere interaction model where the extent of 
the magnetosphere is determined by balancing the outward 
diffusion and inward advection of the stellar magnetic 
field at the inner edge of the disk. They showed that the 
disk-magnetosphere system has two stable torque states 
for certain combinations of the magnetic Prandtl numbers 
and the fastness parameter. If the star is initially spinning 
up, in the absence of extraneous perturbations, the spin-up 
equilibrium eventually vanishes and the star subsequently 
spins down. In its current form, the model does not exhibit 
repeated torque reversals observed. 



The general idea of the bi-state model is that there may 
exist two stable states of the accretion disk correspond- 
ing to the spin- up and spin-down of the NS. If the sys- 
tem is triggered to jump from one state to the other, it 
appears as the torque reversal. How to establish the two 
states is the essential problem for thi s kind of model. In the 
model suggested by lYi et al.l (|l997t ) , the torque reversals 
are caused by alternation between a Keplerian, thin disk 
and a sub-Keplerian, advection-dominated accretion flow 



3. QPOs in 4U 1626—67 and constraints on 
previous models 



The mHz QPOs 

tected with Ginqa (jShinoda et ah 
(lAngelin i et al.'. '1995D, B ^poSA X doi 
RXTE (Komm ers et all Il998t ICha 



4U 1626—67 have been de- 
(iShinoda et all [1991 . ASCA 



wens et al.', '1997 



RXTE 

and XMM- Newton 
iKaur et al.l 



19981 IChakrabartv , ,199i 

Krauss et al. 1, I2007D . More recently, 
investigated the evolution of the QPO 
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frequency in 4U 1626—67 over a long period. It was shown 
that the QPO frequency in 4U 1626-67 during the last 22 
years evolved from a positive to a negative trend: in the 
earlier spin-up era, the QPO central frequency increased 
from ~ 36 mHz inl983 to ^ 49 mHz in 1993, while in the 
subsequent spin-down era, it gradually decreased at a rate 
~ (0.2 ± 0.05) mHzyr"^. However, the lack of observations 
around 1990 does not allow to define an exact time when 
the evolutionary trend of the QPO frequency changed. It 
seems to be somewhat coincident with the torque reversal. 

In accretion-powered X-ray pulsars, the QPO frequency 
is usually regarded to be directly related to the inner radius 
To of the accretion disk. The widely adopted QPO theories 
are the Keplerian frequency model (KFM) and the BFM, 
which consider the QPO frequency as the Keplerian fre- 
quency vk at To and the beat between i/k at rg and the spin 
frequency i's{= ^Is/^tt) of the NS, respectively. In the case 
of 4U 1626-67, we adopt BFM rather KFM as the proper 
interpretation, since the spin frequency is larger than the 
QPO frequency, which means the system would be in the 
propeller phase if we consider the QPO frequency to be the 
Keplerian frequency at rg . 

According to BFM, the frequency at tq is 



and the fastness parameter is 



OJs = 



(10) 



(11) 



If the BFM correctly explains the QPOs in 4U 1626-67, 
we can infer that vk increased during the spin-up era and 
then slightly decreased during the spin-down era. Since the 
the X-ray luminosity Lx originating from disk accretion 
positively correlates with j^k, one would expect a similar 
evolutionary trend in Lx- This is opposite the fact that 
X-ray luminosity of 4U 1 626—67 has been decreasing since 
1977 (jKrauss et al.l . [20071 ). One intermediate implication is 
that the observed change in the X-ray luminosity is not 
due to a change of ma ss accretion rate by the same fac- 
tor (iKaur et al.l . I2008D . For example, the accretion flow 
in 4U1626— 67 might compose multi-components (e.g. disk 
and coronal streams), and disk accretion accounts for only 
a fraction of the X-ray luminosity, but influences the evo- 
lution of the QPO frequency. Alternatively, the X-ray lu- 
minosity variation could be due to a change in obscuration 
by an aperiodically prec essing warped accretio n disk, as 
suggested for Cen X-3 (Raic hur fc Paull . I2008D . A hue of 
evidence in favor of the latter suggestion is that, since the 
discovery of this X-ray source over two decades ago, the 
optical flux from the accretion disk has remained essen- 
tially constant thoug h its overall X-ray flux has declined 
(|Chakrabartvl . ll99a ). 

Now we use the above results to constrain current accre- 
tion torque models proposed for the torque reversals. The 
magnetically threaded disk model always expects that 
gradually increases when the NS evolves from spin-up to 
spin-down, opposite to what we learnt from t he QPO ob- 
serva tions. In the modified-prop eller model (iPerna et al.l . 
|2006[ ) and the bi-state model of iLocsei fc MelatosI (12001 . 
the inner radius of the disk is expected to change between 
a relatively small and large value, corresponding to spin-up 
and spin-down, respectively. If the QPOs are related to the 
behavior in the inner edge of the accretion disk, a transition 



of the QPO frequency is expected. In the model of iYi et all 
(I1997D . the inner radius of the disk might not change, but 
the dynamical change in the disk causes the angular veloc- 
ity at the inner edge to change from uk to Aiy^. Similarly, 
the retrograde disk model suggest that the QPO frequencies 
should vary between uk — Vs and i^k -t- i^s- AH these models 
predict a jump of j^qpo before and after the torque rever- 
sal, which seems to be in contradiction with observations, 
if BFM really works for the QPOs in 4U1626-67. 

4. An alternative explanation 

The above arguments suggest that the torque exerted by 
the accretion disk itself may not play the sole role in ac- 
counting for the torque reversals, and an extra spin-down 
mechanism seems to be required. The latter needs to have 
the following properties: (1) it occurs without requiring 
large decrease in the mass accretion rate (i.e., not the 
propeller-driven AM loss for rapid rotators); (2) it may last 
for long episodes (up to several years); (3) its spin-down 
torque increases with the mass accretion rate. A potential 
candidate is the stellar and disk winds (or outflows) which 
have not been considered seriously in previous works. 

A wind-driven spin-down mechanism for disk-accreting 
magnetic stars has been widely adopted for the classical 
T Tauri stars (CTTS). A large fraction of CTTS are ob- 
served to rotate at approximately 10% of the break-up 
speed, although they have been actively accreting mate- 
rial froin_y3Bir_simoun^^ disks for ~ 10^ — 10^ 
yr (jBertoutl . Il989f ). iKonigll ()1991[ ) applied the GL model 
to CTTSs to explain their slow rotation. However, it was 
pointed out that, when the differentially twisting angle 
between the star and the disk monotonically increases, 
the torque exerted by the field lines first reaches a max- 
imum value, then decreases. This occurs because the az- 
imuthal twisting of the dipole field lines generates an az- 
imuthal component to the field, and the magnetic pressure 
associated with this component acts to inflate the field, 
causally disconnecting the star and the disk (e.g., Aly, 1981; 
lUzdenskv et"all . l2002li . Thus the size of the disk region that 
is magnetically connected to the star is smaller, and the 
magnetic spin-down torque on t he star is significantly les s 
th an in the o riginal GL model (|Matt fc Pudritd . l2005aD . 
M att fc P udritz ( 20 05b. ) further explored the idea of pow- 
erful stellar winds as a solution to the AM problem, and 
showed that stellar winds are capable of carrying off the 
accreted AM, provided that the ratio of the outflow rate 
and the accretion rate is 0.1. In this model a significant 
part of the disk matter is launched as the stellar winds, 
although the mechanism is unknown. If it is due to the im- 
pact of plasma on the stellar s urface from magnetospheric 
accretion streams, the work bv lCranmed ([2008, 2009) sug- 
gested that it could produce T Tauri-like mass-loss rates of 
at least 0.01 times the accretion rate. 

Outflows from the disk-magnetosphere boundary were 
investigated by many authors, both theoretically and nu- 
merically, and episodes of field infiation and outflows were 
obser ved fe.g.lGoodson et al.l . Il997l llQQQt lUstvueova et al.L 
[20061 : IR omanova et al.l . l2009tl . The maximum velocities in 



the outflows are usually of the order of the Keplerian veloc- 
ity of the inner region of the disk. This favors the models 
where the outflows originate from the inner region of the 
disk or from the dis k-magnetosphere boundary. In the most 
recent simulations, iRomanova et al.l (|2009l ) reported that. 
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in the case of slowly rotating stars, the magnetic flux of the 
star can be bunched up by the disk into an X-type configu- 
ration, leading to a conical wincQ, when the turbulent mag- 
netic Prandtl number (the ratio of viscosity to diffusivity) 
> 1 and when the viscosity is sufficiently high, a > 0.03. 
The amount of matter flowing into the conical wind was 
found to be 10 — 30% of the disk accretion rate. 

Observationally there also exists possible evidence for 
the existe nce of winds f rom 4 U1626— 67. With Chandra ob- 
servation ISchulz et al.l ()2002D resolved the Ne/0 emission 
line complex near 1 keV into Doppler pairs of broadened 
{'^ 2500 kms~^ FWHM) lines from highly ionized Ne and 
O, and suggested that they might originate in a disk wind 
driven from the pulsar's magnetopause. The wind mass loss 
rate ~ lO~^°M0yr~^ was shown to be of the same or- 
der as the observed mass accretion rate onto the NS. The 
structure of the emission lines and the helium-like Ne IX 
and O VII triplets support the hypothesis that they are 
formed in the high-d ensity environment of an accretion disk 
(jKrauss et al.l . l2007D . 

In an evolutionary view, perhaps winds (or outflows) 
are inevitably required for 4U1626— 67. It is related to the 
puzzle that the observationally inferred mass transfer rate 
of M ~ 2 X 10^^° MQyr^^ is much larger than theoretical 

expectations M ^ 3 x 10"^"'^ MQyr"-'^ for mass transfer 
from a 0.02Mq donor in a 42 mi nute binary dr i ven b y 
AM loss via gravitational radiation (jChakrabartv I . Il998l) . 
The discrepancy between the measured and predicted mass 
transfer rates indicates that there must be other driving 
mechanisms besides gravitational radiation, and wind mass 
loss is one of the most suitable choices. 

We now propose a model for the torque reversals based 
on the wind/outflow-assisted spin-down mechanism. When 
there is no (or weak) wind from the star and the disk, the 
total torque exerted on the NS comes from the accreting 
material and the magnetic field-disk interaction. The latter, 
however, may not contribu te significantly to the to rque for 
slow rotators, according to lMatt fc Pudrit3 (jlOOS^), so the 
NS experiences a spin-up torque. 



(12) 



During the spin-up stage, the increase of the mass trans- 
fer rate (observed as the increase of the QPO frequency) 
can lead to the bunching of the field lines if the inward 
flow is faster than ou tward diffusion of the field lines 
(|Romanova et al.l . |2009[ ). When the field topology around 
the magnetosphere becomes open, strong stellar -I- disk 
winds are launched, and the NS enters the spin-down stage. 
The total torque becomes 



Nsd = iVo + iVdw + TVs, 



where TVdw and Ns^k are the torques from the disk 
wind and the st e llar w ind, respectively. Calculations by 
iRomanova et al.l (j2009[ ) showed that close to the inner 
boundary of the accretion disk a significant part of the AM 
within the disk matter is carried away by the conical wind. 
Here we assume A^dw = — X-^o with x lying between and 
1. The torque from the stellar winds is 



(14) 



^ In their simulations Romanova et al. (2009) did not take into 
account possible stellar wind. 



where ta is the Alfven radius, at which the poloidal wind 
velocity equals the poloidal Alfven speed. The dimension- 
less factor K takes into account the geometry of the wind 
and is order of unity (k = 2/3 for a spherically symmet- 
ric wind). The magnitude of ta depends on the the mag- 
netic field strength and geometry, mass loss rate, and wind 
sp eeds. A semi-a nalytic, fitted formulation was suggested 
by iMatt &: Pudrit z (2008.) from two-dimensional (axisym- 
metric) MHD simulations. 



R 



(15) 



where K ~ 2.11, m ~ 0.223, R is the stellar radius, B the 
surface magnetic field strength, and Vesc = {2GM / R)^/^ 
the escape speed from the stellar surface, respectively. If 
the wind is launched at the boundary of the magnetosphere, 
then R should be replaced by the magnetospheric radius. 
Adopt typical values of the parameters of 4U1626— 67, i.e., 
M = IAMq, i? = 10^ cm, B = 3 X lO^^ G, ~ O.IM ~ 
^gis gs~^, one can obtain j-q ~ 5.5 x 10 x 10^ cm, and 
fA/fo ~ 1.5 — 6. Assume x ~ 0-5, the ratio of the spin- 
down and spin-up torques can be estimated to be 



I — I 



1.35(^)(^^)(r^ 
^0.5-^^ 0.1 6 



OJs 



)'( — ) 



0.5, 



(16) 



(13) I 



and we get Nsu ^ ~Nsd ^ Ix 10'^^ gcm^s^^, corresponding 
to a spin-up/down rate of ^ (±)8 x lO^^'^ Hzs~^, com- 
patible with the observed values. Note that in the current 
model, the mass transfer rate during the spin-down stage 
can maintain a higher value compared to that during the 
spin-up stage, as shown by the QPO frequency evolution. 
Additionally, Eqs. (13) and (14) indicate that the spin-down 
rate can increases with mass transfer rate, if a roughly con- 
stant fraction of the transferred mass goes into the winds. 

5. Discussion and conclusions 



Based on the measurements of QPOs in 4U1626— 67 and 
the beat frequency interpretation, we proposed a model for 
the torque reversals in this source. The essential idea is 
that the spin-down is induced by stellar and disk winds 
(or outflows) that take away the AM of the NS. Thus a 
significant decrease of the mass transfer rate and possible 
propeller effect are not required. Because of wind mass loss, 
the accretion rate of the NS is not simply the mass transfer 
rate through the accretion disk. Since the latter determines 
the QPO frequencies, it is not expected that there exists 
straightforward, positive correlation between the QPO fre- 
uency and the X-r ay luminosity (or the spin changing rate) 
Bozzo et alll2009D . 

The model seems to be in line with observations of 
other X-ray pulsars besides 4U 1626—67. For example, 
4U 1907-1-09 was found to switch from spin-do wn to spin- 
up without considera bly change in luminosit y ([Fritz et al.L 
2006). Furthermore, lin't Zand et"all ()1998[ ) pointed out 
that, during the spin-dwon stage, the magnetospheric ra- 
dius from the cyclotron line measurements is tq ^ 2400 
km, less than the corotation radius ^ 12000 km, which 
rules out the possibility of the propeller effect being the 
spin-down mechanism. Continuous monitoring of Her X— 1 
showed that its pulse period evolution resembles a saw- 
tooth composed of spin-up and spin-down episodes, and 
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there occurred extremely large spin-down torques up to 
5 times as strong as the spin- up ones, which are very 
likely related to episo dic ejection of matter in Her X— 1 
(jKlochkov et all . I2009D . Signatures of ouflo wing gas was 
also found in the UV spe ctrum of this source (jVrtilek et all . 
[200lt iBorosOTTall . [200l . 

The main limitation of our model is the mechanism for 
the occurrence of the winds, which need s to be explored 
in more detail. In iRomanova et all (|2009f ). strong outflows 
from the inner region of the disk are expected to result from 
field bunching when the mass transfer rate is enhanced, 
and when the magnetic Prandtl number of the turbulence 
is larger than unity. The condition for the bunching of the 
field lines is however, not well understood. It requires that 
the speed of the inward flow of matter in the disk should 
be higher than the speed of outward diffusion of the stel- 
lar field lines. If the accretion rate is determined not only 
by the viscosity but by any other mechanisms of outward 
AM transport, such as by the spiral waves, then this con- 
dition will be satisfied and the bunching of field lines is 
expected. The stellar winds might be driven by some frac- 
tion of the accretion power in the way of accretion shocks 
and/ o r magnetic reconnection events (e.g. lMatt &: Pudr"it3 . 
I2005bl: ICranmeil I2008D . We also note that there could be 
outflows caused by heating of hard X-ray emission of the 
NS, although more likely in the case of spherical accretion 
(jlllarionov fc Kompaneetsl . 119901 ) . These authors suggested 
that, if the X-ray luminosity falls in the region of ~ 2 x 10'^'* 
ergs~^ < Lx <^ 3 x lO'^^ ergs"-^, Compton scattering heats 
the accreted matter anisotropically, and some of the heated 
matter with a low density can flow up and form outflows 
to take the AM away. It is interesting to note that the lu- 
minosity of 4U 1626—67 at torque reversals just fulfills the 
criteria if its distance is ~ 5 kpc. 

The occurrence and disappearance of the stellar and 
disk winds may be accompanied with possible change in 
the structure of the magnetosphere and the accretion disk, 
leading to variation in the radiation features of the NS and 
the disk. Observations with Chandra and XMM-Newton 
showed that the pulse profile of 4U 1626—67 has changed 
significantly from what was found prior to the torque re- 
versal in 1990, sug gesting a change in the geometry of the 
accretion column ( Krauss et al.l . [20071 ) . The X-ray contin- 



uum spectrum was also shown to be closely correlated with 
the torque state. During the 1977—1990 spin- up phase, the 
spectrum was well described by an absorbed blackbody , 
a power law an d a high energy cutoff (|Pravdo et all . Il979t 
iKii et aLl . ll986l) . After the torque reversal in 1990, the time- 
averaged X-ray spec t rum was f ound t o be relatively harde r 



(jOwens et all . Il997t lYi et all . 119971: iKrauss et all . l2007l) 



which is often regarded as an indication of outflows. The 
spectrum was also found to be harder during the new torque 
transition in 2008 than before or after. These results imply 
that the torque reversal is not a simple case of change in 
the mass accretion rate, but there is also a change in the 
accretion geometry in the vicinity of the NS. 

When the inner disk comes closer to the star, there is 
a higher difference in angular velocity between inner disk 
and magnetosphere, and inflation of the field lines is more 
efficient. The difference between the angular velocities can 
lead to cyclic evolution of the field lines - development of the 
toroidal field component, field line opening, and reconnec- 
tion, which were suggested to be accompani ed with energy 
release or flaring activities (e.g. lAlvl . Il985l) . 4U 1626—67 



was indeed seen to fiare dramatically in both X-ray and 
optical on times cales of ~ 100 s before 1990 (Joss et all 
1978; McCli ntock et al.l . 119801: iLi et al.l . [l980). However, 
the re were no fla r ing ev ents seen in any of the observations 
by iKrauss et al.l ()2007t ). The cessation of flaring activity 
may have occurred at the same time as the torque reversal. 
As the NS has entered the spin-up phase since 2008, it is 
interesting to see whether flaring will appear again. 
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